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INCOMPRESSIBLE LIMIT FOR COMPRESSIBLE FLUIDS WITH 

STOCHASTIC FORCING 


DOMINIC BREIT, EDUARD FEIREISL, AND MARTINA HOFMANOVA 


Abstract. We study the asymptotic behavior of the isentropic Navier-Stokes sys¬ 
tem driven by a multiplicative stochastic forcing in the compressible regime, where 
the Mach number approaches zero. Our approach is based on the recently developed 
concept of weak martingale solution to the primitive system, uniform bounds derived 
from a stochastic analogue of the modulated energy inequality, and careful analysis 
of acoustic waves. A stochastic incompressible Navier-Stokes system is identified as 
the limit problem. 


1. Introduction 

Singular limit processes bridge the gap between fluid motion considered in different 
geometries, times scales, and/or under different constitutive relations as the case may 
be. In their pioneering paper, Klainerman and Majda [16] proposed a general approach 
to these problems in the context of hyperbolic conservation laws, in particular, they 
examine the passage from compressible to incompressible fluid flow motion via the low 
Mach number limit. As the problems are typically non-linear, the method applies in 
general only on short time intervals on which regular solutions are known to exist. A 
qualitatively new way, at least in the framework of viscous fluids, has been open by 
the mathematical theory of weak solutions developed by P.-L. Lions [20]. In a series of 
papers. Lions and Masmoudi [21], [22] (see also Desjardins, Grenier [9], Desjardins et 
al. [10]) studied various singular limits for the barotropic Navier-Stokes system, among 
which the incompressible (low Mach number) limit. The incompressible limit is char¬ 
acterized with a large speed of the acoustic waves becoming infinite in the asymptotic 
regime. Accordingly, the fluid density approaches a constant and the velocity solenoidal. 
The limit behavior is described by the standard incompressible Navier -Stokes system. 

In the present paper, we study the compressible-incompressible scenario in the context 
of stochastically driven fluids. Specifically, we consider the Navier-Stokes system for an 
isentropic compressible viscous fluid driven by a multiplicative stochastic forcing and 
study the asymptotic behavior of solutions in the low Mach number regime. To avoid 
the well known difficulties due to the presence of a boundary layer in the case of no-slip 
boundary conditiones (cf. Desjardins et al. [10]), we restrict ourselves to the motion 
in the “flat” A-dimensional torus = ([0,27r]j{o,27r})^: A = 2,3 and on a finite 
time interval (0,T); we set Q = (0,T) x T^. We study the limit as e —>■ 0 in the 
following system which governs the time evolution of the density p and the velocity u of 
a compressible viscous fluid: 

(1.1a) dp -I- div(pu)dt = 0, 

(1.1b) d(pu) -I- [div(pu G u) — uAu — (A -f u)V divu -f ^Vp{Q)]dt = pu) dW. 
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Here p{q) is the pressure which is supposed to follow the 7 -law, i.e. p{g) = where 
7 > N/2\ the viscosity coefficients u, A satisfy 

2 

u > 0, A -I- -u > 0. 

o 

The driving process VF is a cylindrical Wiener process defined on some probability space 
(n, ^, P) and the coefficient ^ is a linear function of momentum pu and a generally 
nonlinear function of density g satisfying suitable growth conditions. The precise de¬ 
scription of the problem setting will be given in the next section. 

The parameter e in (1.1b) is proportional to the Mach number (the ratio of the char¬ 
acteristic flow velocity and the speed of sound). From a physical point of view the 
fluid should behave (asymptotically) like an incompressible one if the density is close 
to a constant, the velocity is small and we look at large time scales. A suitable scaling 
of the Navier-Stokes system results in (1.1b) with a small parameter e, see Klein et 
al. [17]. In the limit of (1. la-1.lb) we recover the stochastic Navier-Stokes system for 
incompressible fluids, that is, 

(1.2a) du - 1 - [ div(u (g) u) — z/Au -|- Vtt] dt = 'f'(u) dW, 

( 1 . 2 b) div(u) = 0 , 

where tt denotes the associated pressure and 'f'(u) = P//<?(l,u), with Vh being the 
Helmholtz projection onto the space of solenoidal vector fields. To be more precise, 
we show that for a given initial law A for (1.1) and the ill-prepared initial data for 
the compressible Navier-Stokes system (1.1), the approximate densities converge to a 
constant whereas the velocities converge in law to a weak martingale solution to the 
incompressible Navie-r-Stokes system (1.2) with the initial law A. This result is then 
strengthened in dimension two where we are able to prove the almost sure convergence 
of the velocities. 

Our approach is based on the concept of finite energy weak martingale solution to 
the compressible Navier-Stokes system (1.1), whose existence was established recently 
in [3] and extends the approach in [12] to the stochastic setting, see Section 2 for more 
details. Similarly to its deterministic counterpart, the low Mach number limit problem 
features two essential difficulties: 

• finding suitable uniform bounds independent of the scaling parameter e; 

• analysis of rapidly oscillating acoustic waves, at least in the case of ill-prepared 
data. 

Here, the necessary uniform bounds follow directly from the associated stochastic ana¬ 
logue of the energy inequality exploiting the basic properties of Ito’s integral, see Section 
3.1. The propagation of acoustic waves is described by a stochastic variant of Lighthill’s 
acoustic analogy: A linear wave equation driven by a stochastic forcing, see Section 3.2. 
The desired estimates are obtained via the deterministic approach, specifically the so- 
called local method proposed by Lions and Masmoudi [21, 22], adapted to the stochastic 
setting. 

A significant difference in comparison to the deterministic situation is the correspond¬ 
ing compactness argument. In general it is not possible to get any compactness in ui as 
no topological structure on the sample space is assumed. To overcome this difficulty, 
it is classical to rather concentrate on compactness of the set of laws of the approxi¬ 
mations and apply the Skorokhod representation theorem. It gives existence of a new 
probability space with a sequence of random variables that have the same laws as the 
original ones and that in addition converge almost surely. However, the Skorokhod rep¬ 
resentation Theorem is restricted to metric spaces but the structure of the compressible 
Navier-Stokes equations naturally leads to weakly converging sequences. On account of 
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this we work with the Jakubowski-Skorokhod Theorem which is valid on a large class of 
topological spaces (including separable Banach spaces with weak topology). In the two- 
dimensional case we gain a stronger convergence result (see Theorem 2.11). This is based 
on the uniqueness for the system (1.2) and a new version of the Gyongy-Krylov charac¬ 
terization of convergence in probability [14] which applies to the setting of quasi-Polish 
spaces (see Proposition A. 4). 

We point out that the gradient part of the velocity converges only weakly to zero 
due to the presence of the acoustic waves, and, consequently, the limit in the stochastic 
forcing $(p, gu)dW can be performed only if is linear with respect to gu. However, 
this setting already covers the particular case of 

$(p, pu) dW = e $1 dW^ + gu $2 dW^ 

with two independent cylindrical Wiener processes and and suitable Hilbert- 
Schmidt operators <l>i and 4>2, which is the main example we have in mind. Here the 
first term describes some external force whereas the second one may be interpreted as a 
friction force of Brinkman’s type, see e.g. Angot et al. [Ij. 

In the case of $(£», gu) = p<I>i, a semi-deterministic approach towards existence for 
(1.1) was developed in [11] (see also [27] for the two-dimensional case). More precisely, 
this particular case of multiplicative noise permits reduction of the problem that can be 
solved pathwise using deterministic arguments only. Nevertheless, it seems that such a 
pathwise approach is not convenient for the incompressible limit. In particular, uncon¬ 
trolled quantities appear in the basic energy estimate and therefore the uniform bounds 
with respect to the parameter e are lost. On the contrary, the stochastic method of the 
present paper heavily depends on the martingale properties of the Ito’s stochastic inte¬ 
gral which gives sufficient control of the expected values of all the necessary quantities. 

We point out that a noise depending on the velocity u in a non-linear way cannot be 
unfortunately handled by the present method. This is due to only weak convergence of 
the velocity due to the oscillations generated by acoustic waves - a problem occurring 
already at the deterministic level, cf. Lions and Masmoudi [21] 

The exposition is organized as follows. In Section 2 we continue with the introductory 
part: we introduce the basic set-up, the concept of solution and state the main results 
in Theorem 2.10 and Theorem 2.11. The remainder of the paper is then devoted to its 
proof. 


2. Mathematical framework and the main result 

Throughout the whole text, the symbols IT*’^ will denote the Sobolov space of func¬ 
tions having distributional derivatives up to order I integrable in . We will also use 
for / £ R to denote the space of distributions v defined on with the finite 

norm 

(2.1) fc^^|cfc(u)|^ < oo, 

keiy 

where Ck denote the Fourier coefficients with respect to the standard trigonometric basis 
{exp(*fcx)}fcgz. 

To begin with, let us set up the precise conditions on the random perturbation of 
the system (1.1). Let (H, (^t)t>o,P) be a stochastic basis with a complete, right- 

continuous filtration. The process W is a cylindrical Wiener process, that is, W{t) = 
^^,> 1 /3fc(t)efc with {/3k)k>i being mutually independent real-valued standard Wiener 
processes relative to {^t)t>o and {ek)k>i a complete orthonormal system in a separable 
Hilbert space it. To give the precise definition of the diffusion coefficient consider 
p £ p> 0, and v £ L^(T^) such that ^/pv £ L^(T^). Denote q = pv and let 
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^(p, q) : U —?> L^(T^) be defined as follows 

<P{p, q)efe = gfc(-, p(-), q(-)) = hfc(u P(-)) + afeq(-): 

where the coefficients afc G R are constants and hfc : x R —?► R. are C^-functions that 


satisfy 


(2.2) 

\ak\^ < oo, 


k>l 

(2.3) 

|hfc(x,p)l2 < C{p^ + \p\' 


k>l 

(2.4) 

^lVphfe(x,p)l2<C(l + lpl^ 


k>l 


Remark that in this setting L^(T^) is the natural space for values of the operator 
pv). Indeed, due to lack of a priori estimates for (1.1) it is not possible to consider 
'P{p, pv) as a mapping with values in a space with higher integrability. This fact brings 
difficulties concerning the definition of the stochastic integral in (1.1) because the space 
L^(T^) does not belong among 2-smooth Banach spaces nor among UMD Banach spaces 
where the theory of stochastic Ito integration is well-established (see e.g. [4], [25], [23]). 
However, since we expect the momentum equation (1.1b) to be satisfied only in the 
sense of distributions anyway, we make use of the embedding L^(T^) 1T“*’^(T^), 

which is true provided / > , and understand the stochastic integral as a process in the 

Hilbert space IT“*’^(T^). To be more precise, it is easy to check that under the above 
assumptions on p and v, the mapping ^(p, pv) belongs to L 2 (H; 1T“*’^(T^)), the space 
of Hilbert-Schmidt operators from il to 1T“*’^(T^). Indeed, due to (2.2) and (2.3) 

llgfc(p>dv)l ]^„,,2 < llgfc(p,pv)l]|i 

fc>l k>l 


(2.5) 


- X! ( / (|hfc(a;,p)l-kplafev]) dx 
<C{p)fN ( r^p-ilhfc(x,p)l2 + ^plafevl2'jdx 


< C'(p)tn [ (p-k p'^-k pjv]^) dx < oo. 


where (p)tn denotes the mean value of p over T'^. Consequently, if 


pGL^nx (0, T), P, dP (g) dt; T^(T^)), 
V^v G L2(H X (0,r),iP,dP® dt;L2(T^)), 


where V denotes the progressively measurable tr-algebra associated to and the 

mean value {p{t))YN (that is constant in t but in general depends on w) is for instance 


essentially bounded then the stochastic integral J^<P{p, pv) dW is a well-defined {^t)- 
martingale taking values in 1T“*’^(T'^). Finally, we define the auxiliary space ilo 3 
via 


iln = <^ u = 


fc>i 


^k^ki 


k>l 


Ikll 


2 

Uo 



k>l 


V — ^ ^ CkGk- 
k>l 


endowed with the norm 
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Note that the embedding it ^ ilo is Hilbert-Schmidt. Moreover, trajectories of W are 
P-a.s. in T];ilo) (see [ 8 ]). 

2.1. The concept of solution and the main result. Existence of the so-called finite 
energy weak martingale solution to the stochastic Navier-Stokes system for compressible 
fluids, in particular (1.1), was recently established in [3]. Let us recall the corresponding 
definition of a solution and the existence result. 


Definition 2.1. A quantity 


IT] 


is called a weak martingale solution to problem (1.la-1.lb) with the initial law A pro¬ 
vided: 

(a) (fl, (.^t), P) is a stochastic basis with a complete right-continuous filtration; 

(b) W is an cylindrical Wiener process; 

(c) the density g satisfies g > 0, t {g(t,-),ip) G (7(0,T] for any if G C°°(T-^) 

P—a.s., the function 1 1 —>■ {g(t, ■),if) is progressively measurable, and 


E 


sup 

tG[0.T] 


< oo for all 1 < p < oo\ 


(d) the velocity field u is adapted, u G L^(fl x (0, T); lE^’^(T^)), 


E 


1*^11 vfl2(t") 


p-\ 


< oo for all 1 < p < oo; 


(e) the momentum gu satisfies t >->• (gu, </>) G C'[0,T] for any (f G C°°(T'^) P—a.s., 
the function 1 1 —> {gu,(j)) is progressively measurable, 

E sup ||gu||^_^ < oo for all 1 < p < oo; 
te[o,T] 

(f) A = Po (e(0),£>u(0))"\ 

(g) for all if G C°°{T^) and ip G C'°°(T^) and all t G [0,r] it holds P-a.s. 

(e’(i),V') = (£'(0),V') + / (gu, Vi/>)ds, 

(giu(t), (p) = (eu(0), ¥?) -I- / (gu^ U,\7(p) ds - V (Vu,Vv 3 )ds 

Jo Jo 

— (A -I- u) [ ( div u, div (p) ds -I- [ (p'^ ,div ip) ds 

Jo £ Jo 

+ f (^(£', £'u)dW,v3), 

Jo 

The present problem requires a refined concept of finite energy weak solution similar 
to that introduced in [2]. The relevant existence result is proved in [3]: 

Theorem 2.2. Assume that for the initial law A there exists M G (0, oo) sueh that 

A|(p,q) G A^(T^) X L^(T^); p > 0, (p)Tiv < M, q(a;) = 0 if p{x) = o| = 1, 

and that for all p G [l,oo) the following moment estimate holds true 


[ ^ 



2 p 


dA(p,q) < Ce- 


LI 










6 


DOMINIC BREIT, EDUARD FEIREISL, AND MARTINA HOFMANOVA 


Then the Navier-Stokes system (1.la-1.lb) possesses at least one weak martingale so¬ 
lution with the initial law A. In addition, the equation of continuity (1.1a) holds also in 
the renormalized sense 


d {b{g), Ip) = {b{g)u, Vip) dt - {{b{g) - b'{g)g) div u, ip) At 

for any test function ip G and any b G c»), b'{g) = for g> gg. Moreover, 

the energy inequality 

p 


E 

sup 

/ 

rl0up , 1 fj-i 'll 

dx 


Lte[o.T] 

JjN 

1 2g J 

. 


( 2 . 6 ) 


+ E 


f [ u|Vup + (A + u)| divup AxAt 
Jo Jt^ 


< c(jp, T) E 
hold for any \ < p < oo, where 


IfN L 


l(e“)(0)P . 1 


da; I + 1 


H{g), H"{g) = 


p'is) 


is the so-called pressure potential. 


Remark 2.3. The pressure potential H is determined up to a linear function. In par¬ 
ticular, one can take 

H{g) = {g'* - 'yg'*~^{g-g) -g'*) 

for any constant p > 0. 

Remark 2.4. Because of (2.6) this solution is called finite energy weak martingale 
solution. The constant c{p,T) depends on T (via GronwalTs lemma) and the constants 
in (2.2)-(2.4) but is independent of e. 


Proof. The proof of (2.6) is somewhat hidden in [3] as multi-layer approximation scheme 
is needed. For the reader’s convenience we give a formal proof on the level of smooth 
solutions. In order to obtain a priori estimates we apply Ito’s formula to the functional 
/(Qj p) = 5 /tn da;. This corresponds exactly to the test with u in the momentum 
equation and in the continuity equation from the deterministic case. We gain 


Z JfN 
_ 1 
“ 2 


p|up da; 


JjN 


l(eu)(0)p 
0 ( 0 ) 


da; — 


0 Jt^ 


I Vup da; da — (A -I- u) 


I divur da; dcr 


10 aT^ 

rt 


A- [ [ pu (g) u ; VudxdtT H—^ [ ( p'*'div udx dtr — - [ [ |u|^dp 
Jo Jt" Jo Jt^ 2 Jq Jjn 

u-^{g,gu)dxdW f g~^ d(^ f pu)dVFy 


JO 

In the following we use the renormalized equation 

rt 


(2.7) 


0 Jr^t 


g^ div u dx do" = — 


/ H{g) 


dx 


/j;N 


H{g{0)) dx 


In fact, (2.7) is a consequence of the mass conservation 


( 2 . 8 ) 



g{t) dx 



p(0) dx 


Vt G [0,T] 
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which holds due (1.1a). Using (2.7) we gain 

1 [ elupdx + uo [ [ |VupdxdCT + ^ [ H(g)dx 

2 JfN Jq JfN Jj3 

1 f f H{g{0))dx 

JfN 


~ 2 


+ 


jN g{0) 

u • $(^i, ^lu) dxdlU + i f g~^ d(^ ( ^(g, gu) dVh^ 


0 

1 
2 


l(pu)( 0 )P 


dx+^ [ H{g{0))dx+ Ti{t)+T 2 {t). 
£ Jt^ 


IfN g(0) 

We apply the g-th power on both sides and then take the expectation. Due to (2.3) and 
(2.8) we have 


(2.9) 


T2(0<^y'/ f g^\gk{g,gu)\'^dxdt<cf f (£»|up + p'>'+ l) dxdt 
2 ^ 7o Jo 2 tn 


< C 


J j ^ + ^i7(e) + l) dxdt. 


As a consequence of Burgholder-Davis-Gundy inequality, (2.3) and (2.8) we gain for 

p > 1 


E 


sup |T 2 (t)| 

p 

= E 

sup 

/ / U • <i>(p, pu) dx dWa- 

-tG(o.r) 


.ieCo,T) 

Jo Jt^ 


= E 

sup 

/ y' / U- g,(£>, £>u)da;d/ 

Jo ^ Jt^^ 


< cE 


< cE 


L^o 


L^O 


^u-g,(g,gu)dxj dt 




IfN 


g |g*(0,£'u)| dx dt 


and by Young’s inequality and a computation similar to (2.9) we gain for every (5 > 0 


E 


1P 


sup |T 2 (t)| 
te(o.T) 


< 5E 


f 

sup / pjupda; 
te(o.T) Jt" 


7((5)E 


_Jo 


+ -^H{g) + 1^ dxdt 


Finally, taking 6 small enough and applying Grownwall’s lemma, the inequality (2.6) 
follows. □ 


Remark 2.5. In Def. 2.1 (j) the continuity equation is stated in the renormalized sense. 
This is part of the existence result in [3] but will not be used in the remainder of the 
paper. 

Concerning the incompressible Navier-Stokes system (1.2), several notions of solution 
are typically considered depending on the space dimension. From the PDF point of view, 
we restrict ourselves to weak solutions (although more can be proved in dimension two), 
i.e. (1.2) is satisfied in the sense of distributions. From the probabilistic point of view, 
we will consider two concepts, namely, pathwise (or strong) solutions and martingale (or 
weak) solutions. In the former one the underlying probability space as well as the driving 
process is fixed in advance while in the latter case these stochastic elements become part 
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of the solution of the problem. Clearly, existence of a pathwise solution is stronger and 
implies existence of a martingale solution. Besides, due to classical Yamada-Watanabe- 
type argument (see e.g. [14], [26]), existence of a pathwise solution follows from existence 
of a martingale solution together with pathwise uniqueness. The difference lies also in 
the way how the initial condition is posed: for pathwise solutions we are given a random 
variable Uq whereas for martingale solutions we can only prescribe an initial law A. 

Note that due to our assumptions on the operator <P, the stochastic perturbations 
that we obtain in the limit system (1.2) is affine linear function of the velocity and takes 
the following form 

^'(v)e/cd/3fe = T’i/I?(l, v)efe d/3fe = (T’ifhfe(l) + afcv)d/3fe. 

Besides, due to (2.2), (2.3) it holds true that 

ll^(v)lli2(U;L2) < + I|v||i2), 

ll^(v) - ^'(w)lli^(ii.i2) < CJlv - wl]|2. 

In dimension three, existence of a strong solution which is closely related to uniqueness 
is one the celebrated Millenium Prize Problems and remains unsolved. Therefore, we 
consider weak martingale solutions, see for instance [7] or [13]. 


Definition 2.6. Let K he a Borel probability measure on Then 


IT) 


is called a weak martingale solution to (1.2) with the initial data A provided 

(a) (n, (.^t), P) is a stochastic basis with a complete right-continuous filtration, 

(b) W is an cylindrical Wiener process, 

(c) the velocity field u is {,^t)-adapted, u G (^^^([O, Tj; L^j^(T^))nL^(0, T; ITJj’.^(T^)) 
P-a.s.^ and 


p 


E 


sup jjuj] 
L(o.r) 


+ E 


Hull 


2 

W1.2(tN) 



< oo for all 1 < p < oo; 


(d) for all ip G C'^„(T'^) and all t G [0,T] it holds P-a.s. 

(u(t), (^) = (u(0), + f (u (g) u, V<p) ds - u / (Vu,Vv3)ds + 


Jo 


Jo 



('f'(u) dW, if). 


Here and hereafter, the substrict div refers to the space of solenoidal (divergenceless) 
functions. 

Under the condition (2.10), the following existence result holds true and can be found 
for instance in [7] and [13]. 


Theorem 2.7. Let K he a Borel probability measure on L‘^{T^) such that for all p G 
[l,oo) 

/ ||vl|^,dA(v)<C(p). 

Jli 

Then there exists a weak martingale solution to (1.2) with initial law P. 


In dimension two, pathwise uniqueness for weak solutions is known under (2.10), we 
refer the reader for instance to [6], [5]. Consequently, we may work with the definition 
of a weak pathwise solution. 


iBy C»([0,T];L2,^(T^)) we denote the space of weakly continuous functions with values in the 
space of divergence-free vector fields in 
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Definition 2.8. Let (.^t),P) be a given stochastic basis with an (^t)-cylindrical 

Wiener process W and let Uq be an -measurable random variable. Then u is called a 
weak pathwise solution to (1.2) with the initial condition Up provided 

(a) the velocity field u is (t^t)-adapted, u G T]; L^j^(T^))nL^(0, T; 

P-a.s. and 


E 


sup 

(O.T) 


1 P 


+ E 


I'^IIwl2(tN) dt 


p-\ 


< oo for all 1 < p < oo] 


(b) u(0) = Up P-a.s., 

(c) for all (f G and all t G [0,T] it holds F-a.s. 

(u(t), (p) = (up, <p)-I- f (u 0 u, Vy) ds - u f (Vu, Vy) ds-I- f {T{u)dW,ip). 

Jo Jo Jo 

Theorem 2.9. Let N = 2. Let (12, (.^t), P) be a given stochastic basis with an 

cylindrical Wiener process W and let Up be an -measurable random variable such that 
Up G LP(fl] T^(T^)) for all p G [1, oo). Then there exists a unique weak pathwise solution 
to (1.2) with the initial condition Up. 

The main results of the present paper are following. 

Theorem 2.10. Let K be a given Borel probability measure on L‘^{T^). Let Ag be a 

2‘7 

Borel probability measure on L'^{T^) x Lt^(T^) such that for some constant M > 0 
(independent of e) it holds true that 




(p,q) e L^(T^^) X L—ir^y, p>—, {p)jN < M, 


p- 1 


<M =1, 


for all p G [1, oo), 


iLjxL 




1^ 

2 p 


dAg{p,q) < C{p), 


LI 


and that the marginal law of Ag corresponding to the second component converges to A 
weakly in the sense of measures on L^^{T^). If Qg, ^e,Wg) is a 

finite energy weak martingale solution to (1.1) with the initial law Ag, e G (0,1), then^ 

Pg—>■ 1 in law on L°°(0, T; L^(T^)), 

Ue-A u in law on (L^(0, T; iT^’^(T^)), w), 
where u is a weak martingale solution to (1.2) with the initial law A. 

Theorem 2.11. Let N = 2 and Up G L^(T^). Let Ag be a Borel probability measure on 
2^7('j'2) ^ (T^) such that for some constant M > 0 (independent of e) it holds true 


that 


13L . 1 


A, (p,q)GL^(T^)xLWT(T^);p> 


M’ 


p- 1 


^ AI, (p)t^ — Af, 


q- Up 


< M)- =1 


If {{n,.^,i.^),F),gg,Ug,W) is a finite energy weak martingale solution to (1.1) with 
the initial law Ag, s G (0,1), then 

Qg^l in L°°{0,T;L'^{T^)) P-a.s., 

Ug —>• u in (L^(0, T; IT^’^(T^)), ir) P-a.s., 

where u is a weak pathwise solution to (1.2) with the initial condition Up. 

Here and in the sequel, the letter C denotes a constant that might change from one 
line to another and that is independent of e. 


^If a topological space X is equipped with the weak topology we write (X, w). 
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3. Proof of Theorem 2.10 

This section is devoted to the study the limit e —>■ 0 in the system (1.1). To this end, 
we recall that it was proved in [3] that for every e G (0,1) there exists 

((^l^^^(^f),P^),Pe,u„lP,) 

which is a weak martingale solution in the sense of Definition 2.1. It was shown in [15] 
that it is enough to consider only one probability space, namely, 

(D^^^P^) = ([0,1],S([0, !]),£) V£G (0,1) 

where C denotes the Lebesgue measure on [0,1]. Moreover, we can assume without 
loss of generality that there exists one common Wiener process W for all e. Indeed, 
one could perform the existence proof from [3] for all the parameters from any chosen 
subsequence e„ at once. The reader is referred to [3] where a similar issue is discussed 
at the beginning of Section 5. 


3.1. Uniform bounds. We start with an a priori estimate which is a modification of 
the energy estimate (2.6) established in [3]. 

Proposition 3.1. Let p G [l,oo). Then the following estimate holds true uniformly in 
e 


E 


sup f (^Q,{t)\u,{t)\^+ ^ {g]{t)-l--f{Q,{t)-l)))dx 

lo<t<TjTN\2 £'^(7-1)" V J 


(3.1) 


+ E 


< UpE 

<Cp. 


£2(7- 1) 

[ [ ulVue]^ + (A + u)l divu, 

Jo Jt^ 

1 




(-Pe(0)|w(0)l" 


+ 


1 


£2(7- 1) 


(6’2'(0) - 1 -7(£’e(0) - l)))da; + 1 


Proof. The first inequality follows directly from (2.6) choosing 
H{g) = (pT' - 'yg'^~^ig - 1) - l) 

and using p(0) < M. Next, we observe that due to the Taylor theorem and our 
assumptions upon A^, it holds 

e/ (p2'(0)-l-7(^?e(0)-l))dx<C£2 

JjN 

and hence the second estimate follows (independently of £). □ 


Corollary 3.2. We have the following uniform bounds, for all p G [l,oo), 

(3.2) Vue G LP{n; L2(o, T; ^^(T^))), 

(3.3) G LP^n-, L“(0, T; ^^(T'^))), 

(3.4) Ue G LP{n-, L2(o, T; 

Proof. The estimates (3.2) and (3.3) follow immediately from Proposition 3.1. For (3.4) 
we use in addition 

[ ge{t)dx= [ £.e(0)da; > ^JT^l 

JjN JjN M 

which is a consequence of (1.1a) □ 
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Let US now introduce the essential and residual component of any function h\ 

h — hggs -j- /ij-es; 

hess = x{Qe)h, X G <7^(0, oo), 0 < X < 1, X = 1 on an open interval containing 1, 

/ires = (1 - x{Qs))h. 

The following lemma will be useful. 

Lemma 3.3. Let P{p) := p'*' — 1 — 7 (p — 1), with p G [0, oo). Then there exist constants 
Cl, C 2 , C 3 , C 4 > 0 such that 

(i) Clip - Ip < P(p) < C' 2 |p- Ip if pG suppx, 

{a) p{p) >€4 if pi supp X, 

{Hi) P{p) > Csp'^ if pi suppx- 

Proof. The first statement follows immediately from the Taylor theorem. The second 
one is a consequence of the fact that P is strictly convex and attains its minimum at 
p = 1. If p ^ suppx and p G [0,1) then the third statement is a consequence of the 
second one. Finally, we observe that the function is increasing for large p S [1, 00 ) 
and its value at p = 1 is zero. This implies the remaining part of {Hi) and the proof is 
complete. □ 


Accordingly, using Lemma 3.3 and (3.1), we obtain the following uniform bounds, for 
all p g[1, 00 ) 

'Qe - ll 


£ J ess 
£>6 - 1 


£ J res 


G LP{n-L°°{G,T-L^{T^))), 
G LP{n-,L°°{Q,T-,L~<{i:^))). 


Therefore, setting := \{Qe — 1), we deduce that uniformly in e 

(3.5) ife e L^’(ll;L“(0,r;L““(T'’2)(T^))). 

As the next step, we want to show that 

(3.6) p, ^1 in LP{^l-L^{0,T-L'^{T^))), 
which in particular leads to 

(3.7) p, G LP{Ll- L“(0,T;L'^(T'^))). 

Then, combining (3.3), (3.4) and (3.7), respectively, we deduce the uniform bounds, for 
all p G [1, 00 ), 


(3.8) p,u, G LP(L!;L“(0,T;(T^^))), 

(3.9) PsUe®u, G LP{n;L^0,T;L^^{r^))). 

Let us now verify (3.6). Since for all <5 > 0 there exists Cs > 0 such that 

p"^ - 1 - 7 (p - 1 ) > Cs\p - IP 


if Ip — 1| > and p > 0, we obtain by (3.1) 




f 

P 

r 

E 

sup 

/ IPe-lpda; 

= E 

sup / l{|g^_i|> 5 }|pe - ipdx 


.0<i<T. 

JjN 




< CiE 


+ E sup / l||g^_i|<5}|pe - Ipda; 

_ 0 <t<TjT^ 

[ (Pe - 1 - XiSe - 1 )) dx 


sup 

0 <t<T Jjn 


+ CS^P < Cse^P + CS'^P. 
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Letting first £ —>■ 0 and then 5 —>■ 0 yields the claim. 

3.2. Acoustic equation. In order to proceed we need the Helmholtz projection Vh 
which projects L^(T^) onto divergence free vector fields 

Moreover, we set Q = Id — Vh- Recall that Vh can be easily defined in terms of the 
Fourier coefficients Ofc (cf. (2.1), in particular it can be shown that both Vh and Q are 
continuous in all IF*’‘^(T^)-spaces, I gM., q € (l,oo). 

Let us now project (1.1b) onto the space of gradient vector fields. Then (1.1) rewrites 
as (using QVf = V/) 

(3.10a) ed(/?e + div Q(£ieU£)df = 0, 

(3.10b) e dQ(peUe) + jVipedt = dt + eQ^{ge, QeU^) dW, 

Fe = vAQUe + (A + u)v div U^ - Q[div(gleUe ® U^)]- \v[q1 - 1 - 7 (£»e - 1)]. 

The system (3.10) may be viewed as a stochastic version of Lighthill’s acoustic analogy 
[18, 19] associated to the compressible Navier-Stokes system. Note that Proposition 3.1 
and (3.9) yield for Z > ^ + 1 using Sobolev’s embedding 

(3.11) F, G LP{0,T;L\0,T;W-^’'^{T^))) 

uniformly in e. 

3.3. Compactness. Let us define the path space X = Xg x X^ x Xg^ x Xw where 

Xg = C^{0,T-L'^{T^)), X^ = {L'^iO,T-W^'\T^)),w), 

= C^([0,T];L^(T'^)), Ah. = C([0, T];ilo)- 

Let us denote by and g,-p(eeUe)^ respectively, the law of Pe, u^, V{ge'n^) on the 

corresponding path space. By we denote the law of W on Ah. and their joint law 
on A is denoted by /i®. 

To proceed, it is necessary to establish tightness of £ G (0,1)}. 

Proposition 3.4. The set {/iuei ^ G (0,1)} is tight on X^- 
Proof. This is a consequence of (3.4). Indeed, for any i? > 0 the set 

BR^{uGL^{0 ,T-,W^’^{T^)y, ||u|U2(o.t;iuu2(t^)) < R} 
is relatively compact in Au and 

1 C 

Mue(-B^) = P(||Ue||L2(o,T;W.C2(TiV)) > R) < ^]E||uE||i2(o,T;M.U2(TiV)) < — 

which yields the claim. □ 

Proposition 3.5. The set {fig^', e G (0,1)} is tight on Xg. 

Proof. Due to (3.8), {div(£i£Ue)} is bounded in L^’(D; L°°(0,T; IA“^’^(T^))) and 
therefore the continuity equation yields the following uniform bound, for all p G [1, oo), 

ge G LP(D;C'°’^([0,r];tA-^’^(T^)). 

Now, the required tightness in follows by a similar reasoning as in Proposition 3.4 to¬ 
gether with (3.7) and the compact embedding (see [24, Corollary B.2]) 

L“(0, T; L'^(T^)) n C°’i([0, Tj; (T^)) 4 ^^([O, Tj; L'>'(T^)). 

□ 
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Proposition 3.6. The set e G (0,1)} is tight on Xg^. 

Proof. We decompose VniBe^e) into two parts, namely, VH{Qe^e){t) = Y^it) + Z'^{t), 
where 

= P^q^(O) - f P//[div(£<£Ue (g) Ue) - i/Aue]ds, 

Jo 

Z^{t)=f VH‘l’igs,gs'Hs)dW(s). 

Jo 

Holder continuity of (Y^). We show that there exists I G N such that for all k G (0,1/2) 
it holds true 


(3.12) E||y®||c»([0,T];W-'.2(TJV)) < C. 

Choose I such that L^(T^) ^ The a priori estimates (3.4) and (3.9) and 


the continuity of P yield 
g 


E 


Y^{t)-Y^{s) 


= E 


^-i,2(TiV) 


f 'P[div{geU^ 0 Ue) 
J s 


e 


^-i,2(TiV) 



/■* 

0 

ft 

< CE 

/ div(£leUe 0 Ue) ds 

+ CE 

/ Au£ ds 


J S 

W-L2(T") 

Js 



ft 

0 

ft 

< CE 

/ PeUe 0 Ue ds 

+ CE 

/ VUeds 


Js 


Js 


uAugjds 
0 

VF-L 2 (T") 


and (3.12) follows by the Kolmogorov continuity criterion. 

Holder continuity of {Z‘^). Next, we show that also 

IE||^®||c''([0,T];W-'.2(T")) < C; 

where I G N was given by the previous step and k G (0,1/2). From the embedding 
L^(T^) W“*’^(T'^), (2.2), (2.3), the a priori estimates and the continuity of Vh we 

get 


E 


Z^{t)-Z-[s) 


= E 


^-i,2(TN) 


PH‘^{ge,ge'J-e)dW 




< CE 


^(&,£'EW)dW 


W- 1 , 2 (JN) 

< CE ( / Ell gfc(£'e,£'£W)||^i dr 


< CE /Ell Sfc (fc; & w) I 


, dr 


fc>i 


fe>i 


< CE 

< 


TiV 


[ge + Peluep + p/)dxdr 


C|t-s| 2 [n-E sup ||y^ Ue||® 2 +E sup < C|t-s |2 

^ 0<t<T 0<t<T ' 

and the Kolmogorov continuity criterion applies. 

Conclusion. Collecting the above results we obtain that 


E||'^H(£'eW)||c''([0,T];W-‘'2(T«) < C" 

for some I G N and all k G (0,1/2). This implies the desired tightness by making use of 
(3.8), continuity of Pff together with the compact embedding (see [24, Corollary B.2]) 

L“(0, T; (T^)) n C'‘([0, T]; ^-'’^(T^)) 4 C„([0, T]; (T^)). 
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□ 

Since also the law/i^u is tight as being Radon measures on the Polish space Xy/ we 
can finally deduce tightness of the joint laws . 

Corollary 3.7. The set {/i®; e £ (0,1)} is tight on X. 

The path space X is not a Polish space and so our compactness argument is based 
on the Jakubowski-Skorokhod representation theorem instead of the classical Skorokhod 
representation theorem, see [15]. To be more precise, passing to a weakly convergent 
subsequence /i® (and denoting by ^ the limit law) we infer the following result. 

Proposition 3.8. There exists a subsequence a probability space (n,.^,P) with X- 
valued Borel measurable random variables (pe, Ug, q^, 114), n S N, and (p, u,q, IT) such 
that 

(a) the law of (g^, Ug, q^, 114 ) given by /i®, e € (0,1), 

(b) the law o/(p, u, q, IT), denoted by /r, is a Radon measure, 

(c) (pg, Ug, qg, ITg) converges P-a.s. to (p, u, q, IT) in the topology of X. 

Let us now fix some notation that will be used in the sequel. We denote by Tt the 
operator of restriction to the interval [0, t] acting on various path spaces. In particular, 
if X stands for one of the path spaces Xg, Xu or Xw and t £ [0, T], we define 

(3.13) rt : X /l[o,t]- 

Clearly, r^ is a continuous mapping. Let (,^f) and (.Tt), respectively, be the P-augmented 
canonical filtration of the process (gg, Ug, iTg) and (p, u, iT), respectively, that is 

= (T((T(rteg, rtUg, rtlTg) U {iV G P(N) = O}), t £ [0,T], 

#t = iT(tT(rtU, rtlT) U {TV G #; P(TV) = O}), t£ [0,r]. 

3.4. Identification of the limit. The aim of this subsection is to identify the limit 
processes given by Proposition 3.8 with a weak martingale solution to (1.2). Namely, 
we prove the following result which in turn verifies Theorem 2.10. 

Theorem 3.9. The process W is a {^t)-cylindrical Wiener process and 

((ii,#,(#4,p),fi,iT) 

is a weak martingale solution to (1.2) with the initial law A. 

The proof proceeds in several steps. First of all, we show that also on the new prob¬ 
ability space the approximations pg, Ug solve the corresponding compressible 

Navier-Stokes system (1.1). 

Proposition 3.10. Let e £ (0,1). The proeess iTg is a {^t)-cylindrical Wiener process 
and 

((fi,#,(#t®),P),Pg,Ug,lTg) 

is a finite energy weak martingale solution to (1.1) with initial law Ag. 

Proof. The first part of the claim follows immediately form the fact that 144 has the 
same law as IT. As a consequence, there exists a collection of mutually independent 
real-valued (.T't)-Wiener processes {f3l)k>i such that iTg = J2k>i Pk^k- 

To show that the continuity equation (1.1a) is satisfied, let us define, for all t £ [0,T] 
and if £ the functional 

L{p,q)t = {p{t),ip) - {p{0),if) - [ (q,VT)ds. 

^0 
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Note that (/o, q) L(p, q)t is continuous on Xg x Xg^. Hence the laws of 
and L(g^, ggUg)t coincide and since (gg, solves (1.1a) we deduce that 

^£^ 5)1 I ~ I =0 

hence {gs,gs%) solves (1.1a). 

To verify the momentum equation (1.1b), we define for all t G [0, T] and ip G C°°{T^) 
the functionals 

Mip,v,q)t = {q{t),p) - {q{0),ip) + [ (q (g) v, V(p) dr - u / (Vv,Vv3)dr 

Jo Jo 

f* a /■* 

— (A + u) / ( div V, div <p) dr H —- / (p'*', div <p) dr 

Jo £ Jo 

Nip,q)t = '^ f (g/c(p,q),<p)^dr, 

k>l 

Nkip,q)t= (g/c(p,q),<p) dr, 

Jo 

let M{p, V, q)s_t denote the increment M{p, v, q)t—M(p, v, q)s and similarly for N{p, q)^^^ 
and Nk{p, q)^,*. We claim that with the above uniform estimates in hand, the mappings 

(P,v,q) H> M(p,v,q)t, (p, v, q) i-S-iV(p, q)*, (p, v, q) H> 7Vfe(p, q)* 

are well-defined and measurable on a subspace of Xg x X^ x Xg^ where the joint law 
of (g, u, pu) is supported, i.e. where all the uniform estimates hold true. Indeed, in the 
case of N(p,q)t we have by (2.2), (2.3) similarly to (2.5) 

f {Sk{p,^),pyds<C'^ f ||gfe(p,q)|||ids<C'. 

M{p,v,q) and Nk{p,v)t can be handled similarly and therefore, the following random 
variables have the same laws 

A^(6’£,Ug,pgUg) - M(pg,Ug,pgUg), 

N{ge, QeUe) ~ iV(p£, feUg), 

J^kiPei P£l^£) ^ ^k{Psi ^£^£)- 
Let us now fix times s,t G [0, T] such that s < t and let 

h ■ 'b’g|[o,s] X ft’u|[o,s] X Xw\[o,s] [0,1] 
be a continuous function. Since 

M(pg,Ug,pgUg)t = / (^>(pg,pgUg)dlT,<p) = ^ / (gfe(pg,pgUg),(p)d/3fe 

Jo 

is a square integrable (.^t)-martingale, we infer that 

[M(pg, Ug, pgUg)]^ - iV(pg, pgUg), M(pg, Ug, pgUg)/3fe - Nkige, PeUg) 

are (^t)-martingales. Besides, it follows from the equality of laws that 

E /i(rgpg, rgUg, TsWe) [M(pg, Ug, pgUg)g,t] 

= E/i(rsPg,rsUg,rgWe)[M(pg,Ug,pgUg)g,t] = 0, 


(3.14) 
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Eh{rsge,rsUi;,rsWe) 


(3.15) 


(3.16) 


[ikf u^, PgUg) t 


= E/i(rs6ie,rsUe,rsll4) 
Eh{rsge,rsUe,rsW^) 


= Eh{rsge,rsUe,rsWe) 

The proof is hereby complete. 


[iH(^£, u^, ^^u^) A^(p£, g^\i^)s i 

[Tf(^£, Ug, A^/j;(^£, g^\i^)s t 


= 0, 


= 0 . 


□ 


Consequently, we recover the result of Proposition 3.1 together with all the uniform 
estimates of the previous subsection. In particular, we find (for a subsequence) that 

(3.17) &^1 in L°°(0,T;T'^(T^)) P-a.s. 

Due to Corollary 3.2 we have the following bounds on the new probability space. 

Corollary 3.11. ITe have the following bounds uniform in e, for all p £ [l,c») and 
I > ^ 

f>e £ LP(fl,L°°(0,T;L™(2,7)(TrW)))^ 

Fe £ LP( 0 ,T;Li( 0 ,r;lT-'’ 2 (T^))) 


where (p^ = ^ ^ and 

Fe = J^AQue + (A + n)V div Ue - Q[div(fcUe ® Ue)]— \v[g2 - 1 - 7 ( 0 ^ - 1)]. 

Proposition 3.12. We have the following convergence P-a.s. 

(3.18) in l 2(0,T;L«(T^)) Vg < ^. 

Proof. Since the joint laws of (gejUi.^Vnigs'^s)) and (p£,U£,q£) coincide, we deduce 
that Qg = Vnige^e) a.s. and consequently it follows from the proof of Proposition 3.6 
that 

(3.19) ]E||7^//(geU£)||c,.QQ^7’].^-i,2(TN)) < C 
for some k £ (0,1) and I £ N. 

Besides, it follows from (3.17) and the convergence of to u that 

(3.20) in L2(0,T;L^(T^)) P-a.s. 

If we pass to the limit in the continuity equation, we see that div u = 0 , which in turn 
identifies q with u. Indeed, due to continuity of P we obtain 

lPff(&fie)^ii in L2 (o^T;L^(T^)) P-a.s. 

Thus with Proposition 3.8 and the compact embedding Lt+t(T'^) 1F“^’^(T^) 

(3.21) VH{geUe)^u in L‘^{0,T;W-^’‘^(T^)) P-a.s. 

Since 

(3.22) divUe^O in L^{0,T; L^(T^)) P-a.s. 


we have also that 

(3.23) Vnue-^u in L^{0,T;W^''^{T^)) P-a.s. 
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Note that (3.22) is a consequence of divu = 0 and the P-a.s. convergence Ug ^ u in 
L^(0, T; c.f. Proposition 3.8. Combining (3.21) with (3.23) we conclude that 

PH{Qe%) ■'PhUe ^ \u\'^ in L^iQ) P-a.s. 


Using Proposition 3.8 yields P-a.s. 


[ {I'PhUeI'^ -Ph{QeUe) -PhUs) dxdt < ||& - l||L~(o,T;L'r)||Ue 
JQ 


2 

L2(0,T;LO 


where s = This implies ||PhUe ||2 —>■ ||u ||2 and hence 

Pud^^d in L‘^{0,T-L^{T^)). 

Combining this with weak convergence in L^(0, T; lU^’^(T^)) (recall Proposition 3.8) 
yields the claim. □ 


In the following we aim to identify the limit in the gradient part of the convective term. 
To this end, we adopt the deterministic approach proposed by Lions and Masmoudi [22]. 
We introduce the dual space 


(T^) = 




In particular, two elements of Wj;*’^(T^) are identical if their difference is a gradient. 


Proposition 3.13. Fori > ^ we have P-a.s. 

div(p£U£ (g) Ue) ^ div(u (g) u) in L^(0, T; Wj7'’^(T^)). 

Proof. Following [22] we decompose 

foUe = U -t Pnigede - u) -t - u) , 

% = u + Ph (ue - u) -I- Q(ue - u) . 

The claim follows once we can show that the following convergences hold true weakly in 


Li(0,r;IU-f(T^))P- 

-a.s.: 

(3.24) 

div (g Ph (ue — u) ^ ^0, 

(3.25) 

div ^u(g) Q(ue - u)^ ^ 0, 

(3.26) 

div (Pnigede - u) 0 ^0, 

(3.27) 

div ^Q(gEUE - u) 0 ^0, 

(3.28) 

div (^Pnigede - u) 0 Ph(ue - u)^ ^ 0, 

(3.29) 

div (^PH{gede - u) 0 Q(ue - u)^ ^0, 

(3.30) 

div ^Q(gEUE - u) 0 Pnids - u)^ ^0, 

(3.31) 

div (^Q(peUe - u) 0 Q(ue - u)^ ^ 0, 


The first four convergences follow from Proposition 3.8, (3.20) and the continuity of 
Ph and Q respectively. The convergences (3.28)-(3.30) are consequences of (3.17) and 
(3.18). In fact, the only critical part is (3.31). First, we need some improved space 
regularity. Similarly to [22] , we use mollification by means of spatial convolution with a 
family of regularizing kernels with a parameter 0 < k <C I. As a matter of fact, thanks 
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to the special geometry of the flat torus T^, the mollified functions can be taken as 
projections to a finite number (which is the smallest natural number > of modes of 
the trigonometric basis {ex.'p(ikx)'\k^z■ In particular, the mollification commutes with 
all spatial derivatives as well as with the projections Vh and Q. For <5 > 0 arbitrary we 
take K = k{5) so small that 


Emn,r - +E||(&iie)" - Q.n.W ^ < S, 

L2(Lt+i) L^(Lf+^) 


(3.32) 

(3.33) E||(&u,)'^-u,"|| 
uniformly in e. We note that the norm E||Ug — Ue|| 


_ 2 j_ +E||u”-Ue|| 2iV 
L2(L-t+i) ® 


lu"^ — u|| 2N < 5, 


2 w can be made uniformly 

small as a consequence of the gradient estimate (3.2). As the mollification commutes 
with div and Q, it is enough to show that P-a.s. 


(3.34) 


div (^Q{g^u^ - Oi") ® Q(fi^ - fi")) ^ 0, 


for fixed k instead of (3.31) (in fact expectation of the (0, T; W)(f^’^(T^))-norm of the 
difference of (3.34) and (3.31) can be estimated in terms of <5 using (3.32) and (3.33)). 
To prove (3.34) we write 

Q(u« - u'^) = - u«) + Q((l - e,-)u^). 

By (3.17), the continuity of Q and the boundedness of u” we know that 

Q((l - e«){i«) ^ 0 in L\Q) 


P-a.s. So (3.34) follows from 

(3.35) div(Q(e?u« - u«) ® Q(p«a« - Q")) - 0, 

in Li(0,r; Wd7'’^(T^)). As div (Qu« (g) Qu'^) = iV|Qu'‘|2, the convergence (3.35) is a 
consequence of 

(3.36) div(Q(&Ue)"® Q(peW)") ^0 in L\0,T-W^I’^{T^)), 

thanks to (3.20) and (3.32). In order to show (3.36) (we need to introduce the function 
div(0eUe) which satisfies Vdie = Q(p£Ue). We have the system of equations 

d{eife) = -V^edt, dV^e = --V(fedt + Fsdt + Q4>(eE, feUE)dM7. 

e 

The right-hand-side only belongs to IF“*’^(T^). So we apply mollification and gain 
4'” = A“^ div(peUe)” and = Q(eEUe)”. The system of equations for <7” and 

reads as 

(3.37) d{eif^) = -A^^dt, dV^^ =dt+ F^ dt + Qd>{g,, g,%)'^dWe. 

We note that for the special choice, where the mollification is taken as the projection 
onto a finite number of Fourier modes, the system (3.37) reduces to a finite number of 
equations. Now, we apply Ito’s formula to the function 


/(£^^,V4'^) = 

with ip G C7)7(T^) arbitrary and gain 

-ifdx 


fjN 




/. 

Jfn 


0 


A4'”V4'g ■ ifdx da — ^ 


0 


(pfV(pf ■ ifdx da 
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+ £ 


0 


■ ipdxAa + e 


/t« Jo 


(pl^ ■ Q$(&, da;. 


And we have 
rt 


0 


■ (pdxdcr 

1 
2 


0 aT" 


V|V^^|2-(pda:dcr- 

0 Jt^ 

= - [ [ : VtpdxdCT, 

Jo Jt’^ 

• V’daidcr = i • 9 ?da;dcr = 0, 

^ Jo Jt’^ 


V^'g (g) V^'” : \7(pdxda 


10 Jt^ 

due to div ip = 0. So we end up with 

V^^(g) : VvJdxdcr =-e / ■ ip dx 

lo Jt^ Jt^ 


+ e 


0 aT-^' 


• 9? da; dcr + £ 


/T"^ ao 


PeV ■ Qe^erdWe dx. 


For fixed n > 0 the right-hand-side vanishes P-a.s. for £ —>■ 0 at least after taking a 
subsequence due to Corollary 3.11, Proposition 3.8 and the properties of the mollification. 
Finally we conclude with (3.36) which implies the last missing convergence (3.31) as 
explained above. □ 

Now, we have all in hand to complete the proof of Theorem 3.9 which implies the 
proof of our main result. Theorem 2.10. 

Proof of Theorem 3.9. The first part of the claim follows immediately from the fact that 
all M4 are cylindrical Wiener processes due to Proposition 3.10. As a consequence, there 
exists a collection of mutually independent real-valued (^t)-Wiener processes {I3k)k>i 
such that W = X]fc>i Pk^k- 

In order to show that (1.2) is satisfied in the sense of Definition 2.6, let us take a 
divergence free test function ip S (7)^^(T'^) and consider the functionals M, N, Nk from 
Proposition 3.10. This way we only study the approximate equation (1.1b) projected 
by Vh and the pressure term drops out. Having (3.14), (3.15) and (3.16) in hand, we 
intend to pass to the limit as £ —>■ 0 and to deduce 


(3.38) 

(3.39) 

(3.40) 


E/i(rsU, rsIT) [M(1 ,u, u)s,t] = 0, 


Eh(rsU, TsW) 
Eh{rsU,rsW) 


[M(l,u, u)^]s.t - iV(l,u)s,t 
[M{l,u,u)Pk]s,t - Nkil,u)s,t 


= 0 , 


= 0 . 


Note that the proof will then be complete. Indeed, (3.38), (3.39) and (3.40) imply that 
the process M(l, u, u) is a (^t)-martingale and its quadratic and cross variations satisfy, 
respectively, 

((M(l,u,u))) = N(l,u), ((M(l,u,u),^fc)) = Nk{l,u), 

and consequently 

|^M(i, a, fi) - ^ (^(1, a) diF, = 0 

hence (1.2a) is satisfied in the sense required by Definition 2.6. 
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Let US now verify (3.38), (3.39) and (3.40). First of all we observe that 
M{ge,Ue,geUe)t ^ M{l,U,u)t a.S. 

due to Proposition 3.8, Proposition 3.13 and (3.21). Application of the Vitali convergence 
theorem together with the uniform estimates (3.4), (3.8) and (3.9) justifies the passage 
to the limit in (3.14) and (3.38) follows. The same argument implies the passage to the 
limit in the part of (3.15) and (3.16) involving M. 

Finally, we comment on the passage to the limit in the terms coming from the sto¬ 
chastic integral, i.e. N and Nk- The convergence in (3.16) being easier, let us only focus 
on (3.15) in detail. As the first step we note that the convergence 

X! {Skigs, &Ue), ^ ® 

k>l k>l 

follows once we show that 


(3.41) (^( 0 e, 0 eUe)-,(p)-)> (^(l,u)-,(p) in L 2 (il;R) P 0 £-a.e. 

To this end, we write 

II (<?(&, fcUe) ■,‘P)- (^(1, U) • , |L,(U;B) 

^k>l ' '^fc>l ' 

= l-l- 

For I 2 we use (2.2) together with (3.21) to obtain I 2 ^ 0 for a.e. (w, t). For Ii we apply 
the Minkowski integral inequality, the mean value theorem, (2.3) and (2.4) to obtain 


k>l 


/i<c( ^l|hfc(&)-hfc(i)||'j <c (^|hfe(fe)-hfc(i)n dx 




k>l 


<C f (l + gj^')|0,-l|da;<c / (l + eJ^Tdx ’’ / 

V / iJjN \ / J L^t 


. JT'v 


\ge - ll'^dx 


where the conjugate exponents p,q G ( 1 , 00 ) are chosen in such a way that 

7 — 1 

p —-— <7 + 1 and 9 < 7 - 
Therefore, using (3.7), (3.17) we deduce 


E / /i dt —7> 0 . 

Jo 

and so for a subsequence / —>• 0 for a.e. (w,t) and (3.41) follows. Besides, since, for all 

P>2, 


E / II (<?>(&, &Ue)-,<p)||^^(,i.R) dr 

J S 

< CE^ 11951122(1 + II95III7 + llv^Uellis) "dr 
^cfl-fE sup IIpeII^t+E sup II < C 

\ 0<t<T 0<t<T J 

due to (3.3), (3.7), we obtain the convergence in (3.15) and therefore u solves (1.2). It 
follows immediately from our construction that for all p S [ 1 , 00 ) 

ueLP(ii;L2(o,r;VF]iV"(T^))). 
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Besides, since we have (due Proposition 3.8 and (3.17)) 

^ u in L^(11;L^((3)) 
lower semi-continnity of the functional 


w I—>■ E 


sup / |w| 

.ie(0,T) 


■dx 


yields u S LP(fi; L°°(0, T; L^(T^))) on account of Corollary 3.11. The usual argument 
about the fractional time derivative (in the distributional sense) implies 

nGLP{n-,C^{[0,nLl,,(T^))) 


and the proof is complete. 


□ 


4. Proof of Theorem 2.11 

In order to complete the proof of Theorem 2.11, we make use of Proposition A.4 which 
is a generalization of the Gyongy-Krylov characterization of convergence in probability 
introduced in [14] adapted to the case of quasi-Polish spaces. It applies to situations 
when pathwise uniqueness and existence of a martingale solution are valid and allows to 
establish existence of a pathwise solution. We recall that in the case of = 2 pathwise 
uniqueness for (1.2) is known (cf. Theorem 2.7). 

We consider the collection of joint laws of 

denoted by /r”’"*. For this purpose we define the extended path space 

= Xg X X^ X Xg^ X Xg X Xu X Xgu X Xw 

As above, denote by the law of W and set u”’™ to be the joint law of 

(Pni ^ Qm 5 ■) 7^ (Pm), IT) on X 

Similarly to Corollary 3.7 the following fact holds true. The proof is nearly identical and 
so will be left to the reader. 

Proposition 4.1. The collection {u"’™; n,m gN} is tight on X^. 

Let us take any subsequence k G N}. By the Jakubowski-Skorokhod theo¬ 

rem, Theorem A. 2, we infer (for a further subsequence but without loss of generality we 
keep the same notation) the existence a probability space (f2,.^,P) with a sequence of 
random variables 

(Pnfc ; 5 Qnfc ; Pm^ : Hm^ j Qm^ : BA;), G M, 

converging almost surely in X"^ to a random variable 

(P, u, q, p,u,4 W) 

and 

^((Pnib I , C(„J. , Pmie ; Umfc , Qmfc ; BA) G ■ ) = U *"(•). 

Observe that in particular, converges weakly to a measure p, dehned by 

p(-) =f’((p,u,q, p,u,q) G ■). 

As the next step, we should recall the technique established in Subsection 3.4. Analo¬ 
gously, it can be applied to both 

(Pnfc j , q»fc,lTfe), (p,u,q,lT) 

and 

(Pmfc ) dmj,, qm,,VF/c), (p,u,4VF) 
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in order to show that (u, W) and (u, W) are weak martingale solutions to (1.2) defined 
on the same stochastic basis where {^t) is the P-augmented canonical 

filtration of (u, u, W^). Besides, we obtain that 

Q = 9 , q = u, q = u P-a.s. 

In order to verify the condition (A.l) from Proposition A.4 we employ the pathwise 
uniqueness result for (1.2) in two dimensions, cf. Theorem 2.7. Indeed, it follows from 
our assumptions on the approximate initial laws A^ that u(0) = u(0) = Uq P-a.s., 
therefore according to Theorem 2.9 the solutions u and u coincide P-a.s. and 

M((e'l,Ui,qi,£>2,U2,q2); (£>l,Ui,qi) = (£l2,U2,q2)) 

= P((e, u, q) = {g, u, q)) = P(u = u) = 1. 

Now, we have all in hand to apply Proposition A. 4. It implies that the original 
sequence (ge, Ug, 7^(^)eUe)) defined on the initial probability space (fl, P) converges in 
probability in the topology of XgX Xg^ to a random variable (p, u, q). Without 

loss of generality, we assume that the convergence is almost sure and again by the 
method from Subsection 3.4 we finally deduce that u is a pathwise weak solution to 
(1.2). Actually, identification of the limit is more straightforward here since in this case 
all the work is done for the initial setting and only one fixed driving Wiener process W 
is considered. The proof of Theorem 2.11 is complete. 

Appendix A. Quasi-Polish spaces 

The so-called quasi-Polish spaces are topological spaces that are not necessarily metriz- 
able but nevertheless they enjoy several important properties of Polish spaces. Let us 
recall their definition introduced in [15]. 

Definition A.l. Let (A, r) be a topological spaee such that there exists a eountable 
family 

{fn ■ X —>■ [—1, Ij; n G N} 
of continuous functions that separate points of X. 

Among the properties of quasi-Polish spaces used in the main body of this paper 
belongs the following Jakubowski-Skorokhod representation theorem, see [15, Theorem 
2 ]- 

Theorem A. 2. Let (A, r) be a quasi-Polish space and let S be the a-field generated 
by {fn] n G N}. If {g,n] n € N} is a tight sequenee of probability measures on {X,S), 
then there exists a subsequence (u/t), a probability spaee with X-valued Borel 

measurable random variables [fk] k G N} and f sueh that /r„j, is the law of fk o-nd 
converges to f in X a.s. Moreover, the law of f is a Radon measure. 

Next, we need to adapt the Gyongy-Krylov characterization of convergence in prob¬ 
ability introduced in [14] to the setting of quasi-Polish spaces. Recall that the original 
argument for the case of Polish spaces follows from the following simple observation 
made in [14, Lemma 1.1]. 

Lemma A.S. Let X be a Polish space equipped with the Borel a-algebra. A sequence of 
X-valued random variables {Y„] n G N} converges in probability if and only if for every 
subsequence of joint laws, {g-nk,mk] k G N}, there exists a further subsequence which 
converges weakly to a probability measure pt such that 

p.{{x,y) € X X X] X = y) = 1. 
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In view of our application in Subsection 4, we are interested in the sufficiency of the 
above condition. 

Proposition A. 4. Let {X,t) be a quasi-Polish space. Let {Yn] n G N} be a sequence 
of X-valued random variables. Assume that for every subsequence of their joint laws 
{/^nfc ,mfc ; k gN} there exists a further subsequence which converges weakly to a probability 
measure fi .such that 

(A.l) qi{{x,y) G X X X; X = y) =1. 

Then there exists a subsequence {Yn,; I G N} which converges a.s. 

Proof. Let / be the one-to-one and continuous mapping defined by 

X {fn{x); n G N}, 

where /„ were given by Definition A.l. Since due to assumption 

in Ax A 

for every (n^), (rufc) and some Y with the law fx, we deduce from the continuous mapping 
theorem that 

{f{Yn,)J{Ym,)) 4 (/(A),/(A)) in [-1,1]" X [-1,1]« 

for every (n^), (rrife). Since [—1,1]^ x [—1,1]^ is a Polish space, Lemma A.S applies to 
the sequence {/(A„); n G N} and the convergence in probability follows. Consequently, 
there exists a subsequence {/(A„,); n G N} which converges a.s. and it only remains to 
prove that A„, converges to A a.s. To this end, we proceed by contradiction: Assume 
that Yni does not converge a.s. to A. Then there exists a set of positive probability 
LI* C LI such that for all uj G LI* there exists a neighborhood Af{uj) of A(w) and for every 
lo G N there exists I > Iq such that A„, (w) ^ J\f{u}). However, as the sequence {/«; n G N} 
separates points of A, there exists n G N such that /„(A„,(w)) ^ /„(A(a;)) and as a 
consequence there exists a neighborhood V(a;) of /(A(a;)) such that /(A„, (w)) ^ V(w). 
This contradicts the a.s. convergence of /(A„,) and completes the proof. □ 
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